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Abstract 



An analytic model of long string network evolution, recently developed by the 
authors, is presented in detail, and modified to describe string loop evolution. 
r^i ■ By treating the average string velocity, as well as the characteristic lengthscale, as 

1) \ dynamical variables, one can include the effects of frictional forces on the evolution 

of the network. This generalized 'one-scale' model provides a quantitative picture 
of the complete evolution of a string network, including the prediction of previously 
unknown transient scaling regimes and a detailed analysis of the evolution of the 
fN| ■ loop population. The evolution of all cosmologically interesting string networks is 

studied in detail, and possible consequences of our results are discussed. 

<N 
(N 

O 1 Introduction 

VD 
Os 

Symmetry breaking phase transitions in the early universe inevitably produce topo- 
logical defects of one form or another. Cosmic strings are of particular interest in 
this context, unlike some other defects, because the evolution of a string network 
does not dramatically alter the standard cosmology. In fact, superheavy strings 
?JZ! ■ associated with a grand unification phase transition provide a much-studied model 
for the initial fluctuations for galaxy formation, also leaving imprints in the cos- 
mic microwave radiation background. But cosmo logical interest in strings is not 
restricted to GUT scales, since strings could have formed at lower energies such as 
electroweak or Peccei-Quinn symmetry breaking with potentially important conse- 
quences, respectively, for baryogenesis or dark matter (axions). Before studying the 
astrophysical consequences of strings, however, one must know how they are formed 
and how they evolve. Due to their statistical nature, the best analytic approach 
consists of doing 'string thermodynamics', that is, describing the string network by 
a small number of averaged quantities. 

The serious analytic study of cosmological string networks began one decade 
ago with Kibble's 'one-scale' model [1] (later modified by Bennett [2]). In this work 
it was assumed that the evolution of the long-string network could be described 
using a single lengthscale, which is usually called the 'correlation length'. One then 
supposes that a scaling solution exists at late times and ends up showing that such a 
solution will in fact exist and be stable subject to conditions on the loop production 
mechanisms. Note that in this model it is conceivable that a string network could 
dominate the energy density of the universe [3]. 
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A step forward in the understanding of string network evolution was provided 
by numerical simulations (see for example [4]). In short, these confirmed the large- 
scale features of Kibble's model, namely regarding the existence and stability of 
the scaling solution (at least in the radiation era), but also showed that it neglects 
important physical processes on small scales. In particular, simulations revealed 
the existence of a significant amount of small-scale structure on long strings, with 
loops being predominantly produced at the smallest scales that can be sampled 
numerically. This change in the understanding of the mechanism of loop formation 
had of course important consequences, notably in the cosmic string scenario for 
galaxy and large-scale structure formation. 

These findings triggered new efforts on the analytical side to try to account for 
small-scale structure. Notably, Austin, Copeland and Kibble developed a model [5] 
where the evolution of the network is described by three different lengthscales, one of 
which aims to explicitly describe the presence of small-scale structure. This model 
also includes a very simple treatment of the effects of gravitational radiation. Apart 
from confirming the predictions of the one-scale model for the large-scale properties 
of the network, the main result of this model is the suggestion that the effects of 
gravitational back-reaction are needed if this 'third lengthscale' is to reach scaling. 
A less attractive aspect of this model is that it has to resort to an unappealingly large 
number of unspecified parameters. There have also been studies of the evolution 
of the linear kink density in what is effectively a 'one-scale model context' which 
anticipated these results by Allen & Caldwell [6] and later by Austin [7]. 

However, it is usually said that it is very difficult to build a house if one starts 
with the roof. It is therefore restrictive to try to build models whose only aim is 
to describe an eventual linear scaling regime. In particular, there is a fundamental 
ingredient in the evolution of a cosmic string network that has been neglected until 
very recently, namely frictional forces due to particle-string scattering, which are 
important for some time after the string-forming phase transition. It should be 
kept in mind that the period immediately after string formation is by no means 
irrelevant, eg for baryogenesis mechanisms involving cosmic strings. Furthermore, 
for electroweak strings the friction-dominated epoch only ends in the matter era, 
lasting almost through to the present day. 

A model of string network evolution including the effects of frictional forces 
has been recently proposed by the authors [8]. Although it should be seen as the 
basis for further work, the model is already predictive enough to be testable in both 
numerical and laboratory experiments, if not cosmologically. 

The model is a simple generalization of the 'one-scale' model in which the 
average string rms velocity becomes a dynamical variable. At present, the model 
does not include small-scale structure (although there are significant hints on how 
to do it) or other potentially important effects such as loop reconnections onto 
the long string network. Nevertheless, it will be shown that this simple model 
provides the first, fully quantitative description of the complete evolution of a string 
network in the early universe (see sections 4-5). In particular, the existence of two 
different transient scaling regimes in the epoch of friction-dominated dynamics are 
established (one of which was previously suggested by Kibble). Also, it is apparent 
in this model that cosmic strings will never dominate the energy density of the 
universe (for reasons other than the statistical physics arguments of Albrecht & 
Turok [9]). 

Cosmic string loops decay fairly quickly after their formation. For this reason, 
their contribution to the seeding of gravitational instabilities or cosmic microwave 
background anisotropies is thought to be subdominant relative to long strings. This 
fact probably explains why the evolution of the loop distribution has been compara- 
tively neglected in the literature. This gap can be filled by the model to be described 
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here. With simple modifications — the more important of them being the use of the 
physical loop size I rather then the correlation length L — this model can also be 
used to study the evolution of the loop distribution. In particular, it will be shown 
that, depending on the parameters characterizing loop production and lifetimes, 
there is more energy density in loops than in long strings. 

Strings (and topological defects in general), of course, are not exclusive to the 
early universe. They exist (and have been seen) in a wide variety of condensed 
matter contexts, including metal crystallization [10], liquid crystals [11,12], super- 
fluid helium-3 [13] and helium-4 [14], and superconductivity [15]. Our generalized 
'one-scale' model can also be used to describe vortex-string evolution in condensed 
matter contexts (with advantages over previously used approaches). In particular, 
some well-known results can be readily reproduced, and new quantitative predic- 
tions regarding loop production can be made. These issues are discussed in a 
companion paper [16]. 

The structure of this paper is as follows. In the next section, after a short 
review of string dynamics, the evolution equations for the 'characteristic lengthscale' 
and the average velocity of the long string network and each individual loop are 
derived and justified. The cases of strings arising from the breaking of gauge and 
global symmetries are both considered. The validity of these 'averaged' evolution 
equations is then tested against simple loop solutions in section 3. Section 4 contains 
a discussion on the importance of the friction force in the early universe, together 
with the analysis of the different scaling laws in the model for both long strings and 
loops. There are also some preliminary comparisons with numerical simulations. We 
then proceed to a detailed and individual analysis of the four physically relevant 
cases: gauge electroweak and GUT, and global axionic and GUT strings are the 
subject of section 5. Finally, section 6 contains conclusions. 

Throughout this paper we will use fundamental units in which h = c = ks = 1. 



2 A generalized 'one-scale' model 



A. String dynamics with friction 

A string sweeps out a two-dimensional surface (the worldsheet) which can be de- 
scribed by spacetime coordinates x M and worldsheet coordinates a a , = x^(a a ) ; 
the line element is then 



where g^ v and 7 a ^ are respectively the four-dimensional spacetime and two-dimen- 
sional string worldsheet metrics. For the case of a gauge (global) string, one can 
then derive the Nambu (Kalb-Ramond) action from the abelian-Higgs (Goldstone) 
model on the assumption that the scale of perturbations along the string is much 
larger than its width 5. (In the global case, one also makes use of the equivalence 
between a real massless scalar field and a two-index antisymmetric tensor field.) 
One finds 



ds 2 = g^ u x^ a x v h da a da b = ^ ab da a da b , 



(2.1) 



S = 




(2.2) 
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where B MM is the antisymmetric tensor field, H^ v \ is its field strength and da^ v is 
the worldsheet area element. Hence the Nambu action is proportional to the area 
swept out by the string. By varying this action one obtains the following equations 
of motion 

, x f 0, Gauge 

-\a a + K x r b x^ a x^ b = ^ H ^ eabxTaX ^ GlQb 8 al . (2.3) 

It should be noted that in the global case \i Q is the 'bare' (unrenormalized) energy 
per unit length. However, it can be shown that if one distinguishes between the 
external and self- field contributions to H and sets H se if = the above equations 
still hold with \i Q replaced by the renormalized energy per unit length, denoted by 
p [17]. 

Still, a crucial ingredient for string evolution is missing. Since strings move 
through a background radiation fluid, their motion is retarded by particle scattering. 
Vilenkin has shown [18] that this effect can be described by a frictional force per 
unit length that can be written 

F f ~ g - " (2.4) 

where v is the string velocity and if will be called the 'friction lengthscale'; its 
explicit value depends on the type of symmetry involved. For a gauge string, the 
main contribution comes from Aharonov-Bohm scattering [19], while in the global 
case it comes from Everett scattering [20]. Then we respectively have 

{^j, Gauge 
^ln 2 (T5), Global ' (2-5) 

where T is the background temperature and (3 is a numerical factor related to the 
number of particle species interacting with the string (strictly speaking, its value is 
slightly different in the two cases, but a common symbol will be used for simplicity). 
Specifically in the gauge case we have 



where this sum is taken over effectively massless degrees of freedom, v a is the phase 
change experienced by a particle transported around the string and b a is 1 for bosons 
and 3/4 for fermions. Hence Aharonov-Bohm scattering will occur for particles with 
non-integer v's; see the paper by Alford & Wilczek in [19] for an example of a model 
with such values. It should also be noted that the Everett scattering formula is only 
valid when the particle wavelength is much larger than the string thickness 5. 

It is then straightforward to show that the frictional force per unit length (2.4) 
can be included in the equations of motion (2.3) by adding the term 

{XT - x^ a Ua) \ { , (2.7) 
(U v being the four-velocity of the background fluid) on its right-hand side. 
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Now consider string motion in an FRW universe with the line element, 



ds 2 =a 2 (r) (rfr 2 -dx 2 ) ; (2.8) 

then U v = (a _1 ,0) and choosing the gauge conditions a = t (ie, identifying con- 
formal and worldsheet times) and x-x' = (ie, imposing that the string velocity be 
orthogonal to the string direction) the string equations of motion with the frictional 
force (2.4) in the background (2.8) can then be expressed as [18,21] 

^('H^-^-Kt)'- (2 - 9) 



e+ 2- + — x 2 e = 0, (2.10) 
• a Ef 



where the 'coordinate energy per unit length' e is defined by 
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and dots and primes respectively denote derivatives with respect to r and a. 1 This 
form of the evolution equations proves to be particularly useful because dissipation 
is naturally incorporated in the decay of the coordinate energy density e, while 
preserving the gauge conditions. 

Incidentally, it has been shown [17] that a global string will behave as a super- 
fluid vortex if it is introduced in a homogeneous background of the form 

tfS = y/pZe ijk (2-12) 

(physically, this corresponds to giving it angular momentum). The interaction 
between this background and the string gives rise to an additional force, known 
as the (relativistic) Magnus force, and (2.9) becomes 

x + (2— + (l — x 2 )x = - f— V + -— xAm, (2.13) 

where 

m = i^x' (2.14) 

Vph 



is the circulation vector; the energy equation (2.10) remains unchanged. 
B. Lengthscale evolution 



We can now proceed to average the string equations of motion to describe the large- 
scale evolution of the string network. We therefore define the total string energy 
and the average rms string velocity to be 




(2.15) 



Note that reparametrizations of <7 can be absorbed into changes of e. 
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Differentiating (2.15) and using (2.10) and (2.16), we see that the total string energy 
density p oc E/a 3 will obey the following evolution equation (in terms of physical 
time t): 

| + (2 H (l + ^) + J)p = 0. (2.17) 

Equation (2.17) incorporates both long strings and small, short-lived loops 
which have (in general) a low probability of interacting with other strings before 
their demise. We shall study the evolution of the long-string network on the assump- 
tion that it can be characterized by a single lengthscale L; this can be interpreted 
as the inter-string distance or the 'correlation length'. Strings larger than L will 
be called long or 'infinite'; otherwise they will be called loops. For Brownian long 
strings, we can define the 'correlation length' L in terms of the network density 2 

Poo as 

Poo = £ ■ (2.18) 

Following Kibble [1] , the rate of loop production from long-string collisions can 
be estimated as follows. Conceptually, we divide the network into a collection of 
segments of length L, each in a volume L 3 . Consider another segment of length 
I moving with a velocity v^; the probability of it encountering one of the other 
segments within a time St is approximately Iv^tjl? . Consistently with our 'one- 
scale' assumption, we then assume that the probability of such an intersection 
creating a loop of length in the range / to / + dl will be given by a scale-invariant 
function w (l/L). The rate of energy loss into loops is then given by 

(dpoo\ Voo f ( £\ £ dl _ _ poo 

where the loop 'chopping' efficiency c is assumed to be constant. Note that in 
previous analyses without friction was assumed to be constant and absorbed 
into the definition of c. 

Finally, by subtracting the loop energy losses (2.19) from (2.17) and then using 
(2.18), we obtain the overall evolution equation for the characteristic lengthscale 
L, 

2^ = 2HL(1 + vl) + ^ + ~c Voo . (2.20) 

Note that with the exception of the expansion term, all terms on the right-hand 
side are velocity-dependent; this will have important consequences (see below). 



C. Loop evolution 



On the other hand, we can also study the evolution of the loop density and distri- 
bution. The traditional approach is to define ne(£,t)d£ to be the number density of 
loops with length in the range (£, £ + d£) at time t; the corresponding loop energy 
density distribution is 

Pi (£,t)d£ = p£n £ (£,t)d£. (2.21) 



Throughout this paper the subscript 'OO' refers to properties of the long ('infinite') string network. 



6 



Note that the total loop energy density is 



Po = J Pl (£,t)d£, (2.22) 

and p = p^ + p 3 From our assumptions on the loop production rate (2.19) it is 
then easy to see that 

£ + ( M(1+o|) + g) w = ,^(<), ( , 23) 

where g is a Lorentz factor accounting for the fact that loops are created with non- 
zero center-of-mass kinetic energy (lost through velocity redshift). However, note 
that this equation is 'static', in the sense that it does not include loop decay mech- 
anisms (eg, via the emission of gravitational, Goldstone boson or electromagnetic 
radiation, as the case may be). 

Instead, we start by using our analytic model to describe the evolution of each 
individual loop. Knowing the energy density transferred from long strings into loops 
and estimating their sizes at formation (see below) , one can numerically determine 
the energy density in loops and other relevant quantities at all times. This formalism 
is that we does not allow for loop reconnections, which are unimportant for GUT 
strings (but can be relevant for high-density electroweak or axionic string networks — 
see section 5); furthermore self-intersections could be included by carefully defining 
an 'effective' loop production size. 

The physical size of a loop is simply given by 



= a(r) f eda- (2.24) 

J loov 



'loop 

its time derivative can be easily calculated using (2.10). However one must still 
subtract energy (hence length) losses due to radiative processes. For the case of a 
gauge string, this can be roughly estimated from the quadrupole formula 

("4) ~g(*°\\ G ,?S, (2,5) 



\dtj rad \dt 3 



(D ~ p£ 3 being the loop's quadrupole moment). Again, note that loop velocity is 
usually assumed constant (v% = 1/2) and not included in (2.25). (This is obviously 
correct in the 'free' regime, but it is not a good assumption in the friction-dominated 
regime.) Then we define 

= -r'Gpv 6 , (2.26) 

0,1 / rad 

where according to numerical estimates T' ~ 8x65 (note that the original parameter 
T was calculated in flat space, where v% = 1/2). Then the evolution equation for 
the physical loop size has the form 

| = (1 - 2vj)m - M _ r'Gpvf . (2.27) 



3 Throughout this paper the subscript 'O' refers to properties of the entire loop population; refers 
to the loops with length in the range {£, £ + d£). 
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Again, all but the expansion term are velocity-dependent. 

For the case of axionic strings, however the emission of gravitational is sub- 
dominant with respect to the emission of axions. The above expressions will still 
be valid with the replacement 

r ' G "^Lrfe< (2 - 28) 

where 5 is the string thickness. 

Now, we will assume that loop production is 'monochromatic', ie that loops 
formed at a time t p have an initial length 

£(t p ) = a(t p ) L(t p ) . (2.29) 

Notice that we are implicitly saying that the loop size at formation depends both on 
the large-scale properties of the network (through the correlation length) and on the 
small-scale structure it contains (through the parameter a). At this stage, since the 
model does not include small-scale structure, we shall resort to physical arguments 
about radiative backreaction etc. to obtain an ansatz for a when necessary. 

With this ansatz the scale-invariant loop production function w (see (2.19)) 
becomes 

w (i ) = - a] , (2.30) 



L J a \L 
and the rate of energy loss into loops becomes 



gixc—r, (2.31) 



Ui / to loops n 

with g being the Lorentz factor as above. 

Hence the energy density converted into loops from time t to t + dt is 

d Po (t)=gpcjjdt; (2.32) 



this corresponds to a fraction 



7— = gc—dt (2.33) 

Poo{t) L 



of the energy density in the form of long strings at time t. Then using our ansatz 
(2.30), the corresponding number of loops produced in a volume V is 

dN(t)=g-^Vdt- (2.34) 
hence the ratio of the energy densities in loops and long strings at time t is 

m m *g __ gL , (t) £ mm^l , 9 ^ t) £ , ( , 35) 

where t c is the moment of the network formation and £(t, t') is the length at time t 
of loops produced at time t' . Here, we neglect the initial loop population at t = t c ; 
this probably has a scale-invariant form (see the discussion in ref. [22]. 
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We can therefore numerically (and, in some simple limit cases, analytically) 
determine the loop density at all times. This generalized 'one-scale' model can 
therefore provide a complete description of a string network. 

D. Velocity evolution 



We now consider the evolution of the average string velocity v. A non-relativistic 
equation can be easily obtained: it is just Newton's law, 

4^-^( 2H + i)- < 236 > 

This merely states that curvature accelerates the strings while damping (both from 
friction and expansion) slows them down. On dimensional grounds, the force per 
unit length due to curvature should be fx over the curvature radius R. The form of 
the damping force can be found similarly. 

A relativistic generalization of the velocity evolution equation (2.36) can be 
obtained more rigorously by differentiating (2.16): 



R V It 



(2.37) 



This is exact up to second-order terms. To obtain the damping term we have taken 
(x ) = (x 2 ) 2 . Writing x 2 = (1 + p • q)/2 (p and p being unit left- and right-movers 
along the string) and defining q = — (p • q) the difference between the two is 

(x 4 )-(x 2 ) 2 = l[(( P .q) 2 )-, 2 ] . (2.38) 

Note that numerical simulations of string evolution indicate that q ra d ~ 0.14 and 
Smat ~ 0.26, so this difference should also be small. As for the curvature term, we 
have introduced R via the definition of the curvature radius vector, 

< 2 - 39 » 

where u is a unit vector and s is the physical length along the string (related to 

1 /2 

the coordinate length a by (is = |x'|<i<7 = (l — x 2 ) eda). The dimensionless 
parameter k is defined by 

((l-x 2 )(x-u)) =kv(l-v 2 ). (2.40) 

The parameter k is related to the presence of small-scale structure on strings: on a 
perfectly smooth string, u and x will be parallel so k = 1 (up to a second-order term 
as above), but this need not be so for a wiggly string. On the other hand, k = for 
a string loop in flat spacetime — in particular, it is easy (although rather tedious) 
to show using its definition (2.40) that this is indeed the case for all Kibble- Turok 
and all Burden loops. In that sense, flat spacetime is the case of maximal small- 
scale structure; this is not surprising, since there is no mechanism available for kink 
decay. In a model including small-scale structure, k would probably be a further 
dynamical variable. For this model, however, we shall use an ansatz interpolating 
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between these two extreme regimes motivated by comparisons with microscopic 
evolution. 

Consider a particular string loop. On large enough scales, they are 'frozen' 
in the background, being conformally stretched by expansion (note that very large 
loops should only form in friction-dominated regimes). Then they should have 
relatively little small-scale structure, and taking k ~ 1 should be a good assumption. 
On the other hand, on small enough scales strings behave as if they were in flat 
spacetime, so one requires that k — > as R — > 0. In particular, in the case of loops, 
demanding that their limiting velocity be vf — 1/2 leads to the requirement that 
k oc R as R — > 0. The remaining point consists in noting that what is dynamically 
meant by 'large' and 'small' scales depends not on the relative size of R and the 
horizon but on the relative size of R and the 'damping length' defined as 



With these requirements in mind, and after comparing with the 'microscopic' (ie, 
unaveraged) evolution of some simple solutions (to be described in the next section) 
one arrives at the following ansatz: 



where x is a numerical coefficient [23] of order unity. Recall that the physical loop 
length is approximately i = 2irR; since we will be considering each individual loop, 
this ansatz applies immediately in that case. 

For the case of the long-string network (for which L = R), the reasoning is 
roughly the same in the regime where L ^> However, the opposite regime never 
arises (at least, in the early universe), so there is no simple way of inferring the k 
behaviour. The closest a network gets is in the linear scaling regime in the early 
universe, where the ratio is less than, but still of order unity. Of course there is also 
the problem that k has a slightly different definition — namely, it is an average over 
the whole network, on a scale of the 'correlation length' L. This makes its physical 
interpretation slightly less clear. For these reasons, when discussing this regime we 
will start by assuming that k is a constant of order unity. Further justification for 
this assumption, and a discussion of the possible use of (2.42) or other alternatives 
can be found in section 4. 

Equations (2.20), (2.27) and (2.37) form the basis of our generalized 'one-scale' 
model, which we will now proceed to apply in several different contexts. We note 
that the velocity-independent 'one-scale' model (2.20) has proved to be success- 
ful in describing the large-scale properties of cosmic string networks in numerical 
simulations. Any deficiencies seem to be associated with the emergence of signif- 
icant small-scale structure, that is, propagating kinks and wiggles on scales well 
below L. In friction-dominated regimes, therefore, we should anticipate improved 
quantitative agreement because of the suppression of this substructure. 

3 'Averaged' versus 'microscopic' evolution 

In order to check the validity of our 'averaged' evolution model, and in particular 
our ansatz for k, we will test it against simple loop solutions. 




(2.42) 
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Firstly, consider a circular loop in flat spacetime but with a constant (non- 
infinite) friction length — ie, a condensed-matter-like situation. We can describe the 
loop trajectory simply by 



x = r(r) (sin cos 0,0) E [0, 2n\ . 
Then equations (2.9,2.10) reduce to 4 



(3.1) 



r+ (1 -r 2 ) 



f 1 



l s r 



+ - =0. 



(3.2) 



k vs. log R/lf 



. 1 



0.01 




Figure 3.1: Two possible ansatzes for the parameter k defined in (2.40), specifically for 
strings in a condensed-matter context as a function of the rescaled loop radius R/£f. For 
the case of string loops, the behaviour at large and small scales has a physical justification 
(see text); numerical matching was improved for circular loops by altering the transition 

point from \ = \/2 (solid curve) to \ ~ 0.57 (dashed curve). However, the ansatz 
corresponding to the dashed curve might be relevant for long strings (see section 4). 



Note that the physical ('invariant') loop radius is R = r/y/l — r 2 , obeying 

R 



R=-r' 



(3.3) 



4 Overdots denote differentiation with respect to conformal time T, rather than physical time t. In 
flat space, T = t. 
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also the 'microscopic' velocity is v = — r and obeys 

^ < 3 - 4 > 

On the other hand, our averaged evolution equations (2.27,2.37) take the form 5 
dR 2 R dv ,„ o, fk(R) v\ , n K . 

Notice the similarity between the two approaches. Loops with size much larger than 
the friction length If will be overdamped, with the velocity being approximately 
given by 

v ~ ^ . (3.6) 




(c) (d) 




Figure 3.2: Comparing the 'microscopic' (solid lines) and 'averaged' (dashed) evolution 
equations for a circular loop in a condensed-matter context. Length and time are in units of 
if, and the time axis is with a logarithmic scale. Plot (a) depicts the log of the (rescaled) 
radius, while (b) depicts the loop velocity. The lower graphs are close-ups of the upper ones. 



In this section, averaged quantities will be denoted by overbars. 
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In this case the two sets of evolution equations actually coincide — hence jus- 
tifying our k = 1 ansatz for large R. As the loop gains velocity r and R become 
significantly different and this equivalence ceases to be valid. When R becomes 
much smaller than if, the loop still looses energy due to friction, but this is no 
longer effective in damping its motion — the loop now begins to oscillate relativisti- 
cally. In particular, over one 'period' v oscillates between and 1 (ignoring nonlinear 
effects near R — due to the finite string width). But we know that the averaged 
velocity should be = 1/2 (in the small-scale limit); this is the physical reason 
why we need k to be a 'phenomenologicaP variable. As we mentioned previously, 
this requirement fixes the behaviour of k on small scales to be as shown in (2.42). 
The remaining question is then how to match the two regimes. 





Figure 3.3: Comparing the 'microscopic' (solid lines) and 'averaged' (dashed) evolution 
equations for a circular loop in a condensed-matter context. Length and time are in units of 
if, and the time axis is with a linear scale. Plot (a) depicts the log of the (rescaled) radius, 
while (b) depicts the loop velocity. The lower graphs are close-ups of the upper ones. 



First of all, we need a clear idea of when (and where) the transition occurs. A 
good guess would be the moment of the 'first collapse', ie, the moment when we 
first have v = 1. In fact, this turns out to be a well-defined event. As was first 
pointed out by Garriga and Sakellariadou [23] (and can be easily seen by analytical 
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or numerical study of the equation of motion (3.2)), circular loops with initial radius 
much larger than the friction length always reach v = 1 for the first time when 



R 



Xc — 0.569 



(3.7) 



col 



Note that £f is the physically relevant case for string dynamics in condensed 
matter contexts (recall that the dynamics in that case is always friction-dominated) . 
Also note that because of friction, all loops will rapidly become (almost) circular. 

After numerically comparing the averaged and microscopic evolution equations, 
we find that the simplest possibility (shown in fig. 3.1), 



k = 



1, 

1 R 

V2if 



f s >x 



(3.8) 



provides excellent agreement for thse circular loops (see figures 3.2-3). In particu- 
lar, this turns out to be significantly better than assuming smoother (and slower) 
transitions between the two regimes. As can be readily seen, this ansatz provides a 
very good fit, considering the lack of parameters available. 




log R/tc 



log t/tc 4 



Figure 3.4: The ansatz for parameter k for gauge GUT string loops, as a function of loop 
radius and time after the formation of the string network, (see Section 5). 

In passing, it is worth pointing out that one can also easily calculate the loop 
lifetime [23]. In the relativistic regime, the R evolution equation can be written 



dR 
~dt 



R 



'2V 



(3.9) 
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so we can immediately estimate that the loop will disappear in a time td ec ~ 2£j 
after its first collapse. 

The case of the circular loop (3.1) in the expanding universe is analogous, with 
the constant friction length being replaced by a time-dependent damping length 
(2.41); also the invariant loop radius is now R = ar/y/i — r 2 . Hence the microscopic 
evolution equations now take the form 

r + (l-r 2 ) + J) =0, (3.10) 



Z-*-K-h)' <3 ' n) 
§- HR ~"t- <3 ' 12) 




Figure 3.5: The ansatz for parameter k for gauge electroweak string loops (compare with 
figure 3.4). Note that these units t eq ~ 22.1, friction ceases to dominate the dynamics at 
t* ~ 25.8 and the present time is t a ~ 28.5. 



Neglecting the gravitational radiation term (whose form has been established else- 
where), the averaged evolution equations are 
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Corresponding to this change, we simply modify if to id in our ansatz for k, 
which becomes (2.42) (see figures 3.4-5). However, the numerical value of the ratio 
of R and id when the loop first collapses, x, is now slightly smaller: 

_ J 0.431 , Radiation , . 

X= \ 0.380, Matter ' ^ ' 

Numerically we find that a slightly larger value, x ~ 0-5 provides the best fit in 
both cases. It is also interesting to note that while the radiation-era evolution is 
rather insensitive to the value of x ( m the range x ~ 0.3 — 0.6 say), the matter-era 
evolution exhibits a stronger dependence. 
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Figure 3.6: Comparing the 'microscopic' (solid) and 'averaged' (dashed) evolution equa- 
tions for the physical radius (a), the radius relative to the horizon (b), and velocity (c) of a 
circular gauge GUT loop formed at t = t c with radius R = 10£ c ; (d) is a close-up of (a). 
Time is in orders of magnitude from the moment of loop formation. Radiative mechanisms 
are not included. 



Figures 3.6-7 depict the evolution of a loop with initial radius i?j = 10£j in the 
radiation (ti = t c ) and matter (ti = t eq ) eras — see section 4. Again, the loops are 
initially overcritically damped, the approximate velocity being 

(3-15) 
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The effect of damping is essentially twofold. Firstly, it delays the moment when the 
loop first collapses. While its velocity is non-relativistic, there is no loss of length 
through velocity redshift (see (3.13)), and so the physical loop radius can grow to 
a size much larger than the initial radius. As it picks up speed, however, it starts 
losing more energy. As we mentioned above, the first collapse still occurs for R ~ id] 
in the relativistic regime, the loop loses energy at each oscillation as before. 

Finally, when friction has switched off (and the period of oscillation is much 
shorter than the expansion rate) the loop starts to oscillate with constant physical 
amplitude (see (3.13)). Notice that due to the additional effect of friction this energy 
loss is much larger in the radiation era. In this case the 'final' physical radius is 
almost equal to the initial radius, whereas in the matter era it can be more than one 
order of magnitude larger. When this final stage is reached, gravitational radiation 
or other preferred decay channels cause the loop to shrink further. 
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Figure 3.7: Comparing the 'microscopic' (solid) and 'averaged' (dashed) evolution equa- 
tions for the physical radius (a) and the radius relative to the horizon (b) of a circular gauge 
GUT loop formed at t = t eq with radius R — 10t eq . Time is in orders of magnitude from 
the moment of loop formation. Radiative mechanisms are not included. 

Having thus established, in simple but physically relevant cases, the validity 
of our generalized 'one-scale model', and in particular of the ansatz for k, we now 
proceed to apply it to the study of cosmic string evolution; the study of string 
evolution in condensed matter contexts is left for a companion paper [16]. 



4 General scaling results 



A. Introduction 



In the early universe the friction lengthscale increases with time, so friction will 
only be important at early times [23]; however, the meaning of 'early' is, as we will 
see, model-dependent. Again let T c be the temperature of the string-forming phase 
transition; the corresponding time of formation is 
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where / is given by 



(4.2) 



and M is the number of effectively massless degrees of freedom in the model (e.g., 
M = 106.75 for a minimal GUT model, but it can be as high as 10 for particular 
extensions of it). Then in the case of a gauge symmetry breaking the friction 
lengthscale can be written 



1 t 



3/2 

3\3/2 i 



Radiation era 



(!) 



X72, Matter era 



(4.3) 



and for the case of a global symmetry 
i^ln(i)[ln(f^)] 2 , 



In 



+ +1/3 
8 t <= t eq 
A t4/3 



Radiation era 
Matter era 



(4.4) 



The constant 6 is a measure of the importance of the friction term in the evolution 
equations; its value is 



6 



tr 



A 
v7 Ktpi 



1/2 



(4.5) 



Note that as we mentioned previously, (3 is not exactly the same in the gauge and 
global cases. Also, (4.3) is valid almost immediately after t c (so we will begin 
studying the evolution of the network at ti t c ), but (4.4) is only valid for t/t c > 
6/ A. The string energy per unit length can again be written 



ix = 



rp2 



Gauge case 



T c 2 ln(f), Global case ' 



(4.6) 



Defining t* as the time at which the two damping terms in (2.9) and (2.10) 
have equal magnitude we find 



166> 2 (lnf)" 



In 



( 6 *c\ 

\\u) 



Gauge case 
Global case 



(4.7) 



provided this is still in the radiation era; otherwise, in the matter era we obtain 

1/2 



(i) 1/2 *(£ 

4 (l) 1/2 ^(^) 1/2 ( in f) _1 



In 



8 / ^eg 

A I t c 



1/3 



4/3 



Gauge case 



Global case 



String dynamics is friction-dominated from t c until t*, after which motion will 
become relativistic or 'free'. A simple heuristic argument (see, for example [24]) 
due to Kibble suggests that in the damped phase the correlation length will scale 
as L oc t 5 / 4 ; it is worth noting that this argument is rather similar to the ones used 
in condensed-matter contexts. 
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We now address the problem of initial conditions. Starting with the correlation 
length, by causality this must obviously be smaller than the horizon; but on the 
other hand, it must be greater than the friction lengthscale since friction is initially 
dominating the dynamics. We will therefore assume (following [24]) that 

l H <U<U. (4.9) 

For simplicity, we now concentrate on the gauge case, and define L = L/t c . Then 
(4.9) can be written 

6- 1 <L i <l; (4.10) 

these two extreme limits could correspond to a rapid second-order phase transition 
(Li ~ 9~ l ) or a slow first-order transition (Li ~ 1). Now, the parameter f3 in (2.6) 
can be written 

£=^1^, (4.11) 

(where < u> < 1) and using the definitions of 9 and /, one finds that the initial 
ratio of the string and background densities obeys 




^6003)^ (< .75). (4.12) 



In fact, by analyzing the evolution equations (2.20,2.37) one can show that for 
all physically reasonable values of c and k these bounds hold for all subsequent times. 
Consequently, cosmic strings can never dominate the universe. It can analogously be 
shown that this is also true in the case of a global symmetry, although the required 
algebra is slightly more complicated. Note that this is an entirely different (and 
independent) argument to Albrecht & Turok's [9] statistical physics argument — in 
a sense, it is a 'thermodynamical' argument. 

As for the initial velocity, it can be estimated by a very simple (and rough) 
argument. The friction lengthscale can be approximately interpreted as the distance 
that a piece of string can travel before it is stopped by the friction force. Then by 
comparing the energy contained in strings in a volume L 3 and the work done by 
the friction force in stopping them we obtain 



t'U~^Ll (4.13) 



and substitution of t{ yields 



As expected, highly-curved strings will have high velocities and conversely strings 
in low-density networks will have small velocities. 



B. Long-string scaling laws 



Analysis of the evolution equations for the correlation length and velocity of the 
long string network (2.20,2.37) reveals the existence of three types of scaling regimes, 
which we now describe in detail. Two of these regimes are transient, occurring in 
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the friction-dominated epoch. In this situation strings should have very little small- 
scale structure, and we should have k = 1. In analogy with our discussion for string 
loops in section 3, we should also find that the string velocity in the radiation era 
is 

^OCy, (4.15) 
±J 

which is in fact the case (see below). The third is the well-known linear scaling 
regime. 

Recall that t c (defined in (4.1)) denotes the time of string formation; on the 
other hand, t s will denote the time at which the relevant period of evolution starts. 

Stretching regime 

This is a transient regime that occurs in the beginning of the friction-dominated 
phase, provided the initial string density is sufficiently low. With t s = t c we get 



1/2 



(4.16) 



v = 



t 

t 



In 



A t_ 

6 t s 
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Gauge case 
Global case 



(4.17) 



The reason why this regime arises is physically obvious. If we start with a 
small string density, the correlation length will be close to the horizon, and much 
larger than the damping length. Hence long strings are fairly straight and have very 
little small-scale structure. This is therefore analogous to the situation in condensed 
matter — the strings have very small velocities and are 'frozen', being conformally 
stretched by the expansion. 

Mathematically, the L oc t 1 ^ 2 law comes immediately from the fact that the 
friction, loop production and redshift terms in the evolution equation for L are all 
velocity-dependent, and can therefore be neglected when compared to expansion. 
Note, however, that unlike the condensed matter case there is no additional loga- 
rithmic correction to the L oc t 1 / 2 law in the case of a global string network. Finally, 
and although this is not cosmologically relevant, it should be pointed out that the 
corresponding regime in the matter era would be 



Loct 2 / 3 , voct 4 / 3 



(4.18) 



Kibble regime 

This is a transient regime that also occurs in the friction-dominated epoch (following 
the stretching regime when such a regime exists). In this case the scaling laws are 



L 



2(1+5) 
361 

2(1+5) 
3d 



1/2 / x 5/4 



1/2 



Gauge case 

(^ 5/4 KH)( ln T) 1/2 > Globalcase 



(4.19) 



26(1+5) 
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20(1+5) 



V2 / t \ 1/4 



1/2 



Gauge case 



(£) V4ln (H)( ln T) 1/2 > Globalcase 



(4.20) 
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This is a high-density regime, arising when the correlation length L is close 
to the friction length — either because it started that way or because it becomes so 
during a period of L oc t 1 / 2 evolution (recall that in the radiation era the friction 
length grows as if oc t 3 / 2 ). 

Although friction still dominates the dynamics, the higher string curvature 
(and consequently higher velocity) means that the network is now chopping off a 
considerable amount of energy into loops — in fact, proportionally more than in the 
final, linear regime. However, note that there is still no small-scale structure (ie, 
we still have k ~ 1). Therefore this shows that small-scale structure is not the 
only determining factor for loop production — the long-string density (ie, L) is just 
as important. This is the reason why we chose the form of our loop production 
ansatz (2.29), explicitly separating the effects of large (through L) and small-scale 
structure (through the parameter a). 

The L oc t 5 / 4 scaling law was previously suggested by Kibble. However, note 
that the above relations only hold in the radiation era. In the corresponding regime 
in the matter era for very light strings, we have the scaling laws 

Loct 3/2 , voct 1/2 . (4.21) 



Linear regime 

This is the well-known 'final' scaling regime — it is always the endpoint of cosmo- 
logical string evolution, arising when the friction lengthscale becomes subdominant 
with respect to the Hubble length. Assuming that k (as well as c) are constant in 
each regime, we find the following scaling laws 



L = 



[k r (k r + Cr)] 1 ^ 2 t, Radiation era 



{krn ~~\~Cm ) 



1/2 



t, Matter era 



(4.22) 



v = 



1/2 



\k r ~\~c r ) 

fern 



1/2 



Radiation era 



Matter era 



(4.23) 



Now, the simplest way to proceed is to look for the values of c and k that match 
the simulations: 

c r ^0.24, fej.w0.18, (4.24) 



c m ^0.17, fe m ^ 0.49. (4.25) 

Hence, according to our interpretation of this model, it predicts a larger loop pro- 
duction rate and more small-scale structure in the radiation era (recall that more 
small-scale structure corresponds to a smaller fe) — which is exactly what is seen 
in numerical simulations [4]. This shows that our interpretation of k as being re- 
lated to the presence of small-scale structure is at least qualitatively correct. Notice 
that the scaling parameters are much less sensitive to variations in c than those of 
previous analytic models; this will be relevant below. 

Furthermore, we can get some feeling for the validity of our ansatz (2.42) for k 
by finding out which values we obtain with the known values of the ratios ( = L/t 
in the radiation and matter eras [4]. In the first case we have ( r ~ 0.28 , k r ~ 0.19 in 
excellent agreement with (4.24). In the matter era the scaling value of £ m ~ 0.535 
also gives good agreement with the simulations, since fe m = 0.50 with % = 1 (the 
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dashed line ansatz in fig. 3.1). Note, however, that matter era correspondence is 
upset if we use the x = 0-57 ansatz modified to match collapsing circular loops. 
Nevertheless, these results indicate that the ansatz (2.42), originally justified for 
string loops, is probably extensible to the long string network, although the point 
X where the transition between the constant and linear regimes takes place should 
be studied in greater detail. If this is indeed the case, then it appears that this 
simple model goes a long way towards solving the 'matching' problem at the matter- 
radiation transition — previously, a significant failing of the 'one-scale' model. 

Of course, these results are a manifestation of the fact that we require addi- 
tional degrees of freedom to incorporate small-scale structure satisfactorily (see, for 
example, ref. [5]), but perhaps there are fewer required than at first anticipated. 

C. Loop scaling laws 

We already know how to determine the (relative) loop density at all times; we 
simply have to evaluate 

^=^jlW) l -W) M ' (426) 

where t c is the network formation time and £(t, t') is the length at time t of loops 
produced at time t' . If necessary, we can also analyze the distributions of loop 
lengths, etc. — although this will not be done in this paper. There is still, however, 
one point to be discussed — we must propose an ansatz for the parameter a, which 
should be related to the presence of string small-scale structure. 

We know on physical grounds that a ~ 1 in the friction dominated epoch, since 
strings are non-relativistic and any wiggles are quickly erased; this conclusion is also 
supported qualitatively by observing the actual evolution of networks in quenched 
liquid crystals [12]. On the other hand, all studies of radiative backreaction lead 
us to expect a to be a constant (significantly less than unity) in the linear scaling 
regime. For example, in the case of a gauge GUT string network numerical simula- 
tions [4] show that olgut < 5 x 10~ 3 (note that our definition of a differs from that 
in the numerical simulation papers — where it is defined as £/t). We will therefore 
use the simple ansatz 

1 + a sc -r 

a(t) = -— 1 ±, (4.27) 

' t* 

where a sc is the constant scaling value. Note that this is such that in the transition 
between the damped and free regimes the physical loop length at formation (£ = aL) 
is constant. In this simple way we can 'phenomeno logically' account for the build- 
up of small-scale structure as seen in the numerical simulations [4]. In particular, 
for the case of GUT strings, this build-up takes about 4 orders of magnitude in 
time, in agreement with a result obtained by Allen & Caldwell and by Austin [6] by 
estimates of the evolution of the linear kink density in a 'one-scale' model context. 

Although (4.26) cannot be evaluated analytically in general due to the com- 
plicated behavior of the integrand (and in particular of the £(t,t') factor), it is 
possible (under some simplifying assumptions) to obtain an analytic solution in the 
linear regime. For simplicity we will neglect the loops existing at the start of this 
regime, t = t s . Since loops are much smaller than the horizon, we can approxi- 
mately neglect the effect of expansion; furthermore, we can also assume that the 
loop velocity is close to the terminal velocity, vg ~ 1/V2, so that we approximate 
T'v 6 ~ T. These two simplifying assumptions mean that our calculation will be a 
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Figure 4.1: The interval in t giving a non-zero loop length contribution to the loop 
density integral (4.26) for different times t. 



slight underestimate of the true loop density. The loop decay function in this case 
is the following 

£(t,t') = [ e - T Gv(t-t'), t'<t<t' + T ^ ^ (428) 

[ 0, Otherwise 

with the initial loop length being 

£' = a(t' , (4.29) 

and we used ( = L/t. Note that the coefficient TGn will be altered if the pre- 
ferred loop decay channel is not gravitational radiation, but Goldstone bosons or 
electromagnetic radiation. With these assumptions, (4.26) reduces to 

= ^MW, (4.30) 

"set, J f(t) 1 

where 

f t a , t < % 

and 
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rho vs . log 



t/ts 




Figure 4.2: Analytical estimate of the evolution of the ratio of the loop and long string 
energy densities in the linear scaling regime, for a = 10 3 and TGfJ, = 65 X 10 6 , in 
the radiation (solid) and matter (dashed) eras. Time is in orders of magnitude from the start 
of the regime. 



The path of integration for all times is schematically represented in figure 4.1. 
The integral can be therefore conveniently split into two cases, the transition point 
corresponding to the moment of the decay of the loops formed at t s . Then one finds 
the following solutions 



Pojt) _ gcvpc 
Poo(t) a sc C 



X 



\(Ca sc + TGn) (fj - l) - \YGn (fj - 
i(C« sc + FGfi) (jz - l) - \TGp - 1 



t < ¥ 



t > 



(4.33) 

Again, this can be confirmed numerically (see section 5). Figure 4.2 depicts a 
typical situation; after an initial build-up while no loops have decayed, we reach a 
constant ratio. Note that (4.33) should be testable against numerical simulations. 
With the values of c and k quoted above and assuming that a = 10 -3 and TG/1 = 
65 x 10~ 6 this model predicts the scaling density ratios to be 



Poit) 
Pooit) 



2.1, Radiation era 
2.3, Matter era 



(4.34) 



Notice that it is possible (at least in principle) to measure both a and the 
ratio of the densities from numerical simulations (r can also be estimated in this 
way) . Hence the prediction of this model for the density ratio in the scaling regime 
(4.33) and even the approach to it can be tested numerically. Figure 4.3 shows 
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log rho vs . log alpha 
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Figure 4.3: Expected ratios of the loop and long string densities in the linear scaling 
regime, in the radiation (solid) and matter (dashed) eras for gauge GUT strings, in terms of 
the small-scale structure parameter Ot and with YGfl = 65 X 10 — 6 . 



the predicted ratios for gauge GUT strings in the radiation and matter eras, at 
fixed TG/i for all possible values of a (the upper limit is inferred from numerical 
simulations, the lower limit assumed on physical grounds). Note that if loops are 
large enough, then their density can be dominant (by a factor of up to 50) and the 
ratio will be even higher in the matter era, converse to the opposite limit for small 
loops. The reason is that in the matter era loops live longer (see section 3), a factor 
compensated by the lower loop production rate. 

Finally, let us just briefly mention two other helpful quantities in characterizing 
the evolution of a cosmic string network. The first is the fraction of the energy 
density in long strings at a time t which is converted into loops within the next 
Hubble time; this is given by 

/H = 4lT' (4 - 35) 

and the limiting values for GUT strings in the radiation and matter eras are respec- 
tively 

. J 0.81, Radiation era , 
}h ~ \ 0.18, Matter era ' ^ ' 

The second is the number of loops chopped off by the long string network per Hubble 
volume per Hubble time, given by 

c v, 
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with asymptotic values 

_ \ 2.2 x 10 5 , Radiation era „o\ 
U ~ \ 4.5 x 10 3 , Matter era K ■ ) 

These quantities should also be measurable in numerical simulations. 



Properties 


Gauge EW 


Global Ax. 


Global GUT 


Gauge GUT 


TjGeV 


10 2 


10 10 


10 15 


10 ie 


G/i 


1Q -34 


10 -18 


10" 8 


10" 6 


A 


Irrelevant 


0.2 


0.2 


Irrelevant 


e 


2 x 10 15 


1 x 10 8 


1 x 10 3 


29 


tc/tpi 


4 x 10 32 


3 x 10 16 


3 x 10 6 


3 x 10 4 


ti/t c 


1 


60 


60 


1 


U/t c 


7 x 10 25 


2 x 10 8 


351 


855 



Table 5.1: Some relevant scales for the evolution of the four cosmologically interesting 
string networks. 



Case 


Gauge EW 


Global Ax. 


Global GUT 


Gauge GUT 


Solid 


10" 1 


0.9 


0.9 


0.9 


Dashed 


10" 4 


5 x 10" 2 


0.2 


0.2 


Dotted 


10" 8 


10" 3 


0.02 


0.05 



Table 5.2: Initial conditions on the ratio of the lengthscale to the horizon, to be 

used in the figures and discussion below. 



5 The cosmologically relevant networks 



A. Introduction 



We can now discuss the four cosmologically interesting cases : gauge electroweak 
and GUT, and global axionic and GUT strings. To summarize the cases to be 
studied numerically, some relevant quantities in each case are listed in table 5.1. 
The initial conditions taken for the correlation length (expressed in terms of the 
ratio L/t) are summarized in table 5.2. Note that the initial condition on velocity 
is not independent — see (4.14) in section 4. 
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(b) 
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Figure 5.1: The complete evolution of an electroweak gauge long-string network. Plots 
successively represent the exponents of the power-law dependence of L (a) and V (b), V 
itself (c), and the ratio of the long string and background densities (d). The horizontal axis 
is labeled in orders of magnitude in time from the moment of string formation, and initial 
conditions are specified in table 5.2. 

There is another parameter that we must estimate, namely the typical loop 
size at formation in the linear scaling regime, a 8C . As we have seen in the above 
discussion, this is crucial in determining the ratio of the energy densities in loops 
and long strings. We will assume that a sc ~ G\i. 

B. Electroweak and axionic strings 

We start by noting that for electroweak and axionic strings we cannot meaning- 
fully treat the case of very high string density (that is, the case where L is initially 
close to the friction lengthscale). This is because in this case loop reconnections 
onto the string network may have a significant effect. In all other situations the 
effect of such reconnections can be neglected. 

For gauge electroweak strings (see figures 5.1-2) the epoch of friction-dominated 
dynamics ends well into the matter era. One of the most important results on this 
paper is that the early-time evolution of the network depends crucially on the initial 
string density. Note that among other things this depends on the order of the phase 
transition at which the strings form. If this density is low, strings have very small 
velocities and according to our previous discussion the network starts evolving in 
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(a) (b) 
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Figure 5.2: Characteristics of loop production by electroweak strings: log plots of the 
fraction of the energy density going into loops per Hubble time (a) and of the number of 
loops produced per Hubble volume per Hubble time (b). The horizontal axis is labeled in 
orders of magnitude in time from the moment of string formation, and initial conditions are 
specified in table 5.2; we have taken a ~ 10 — 31 and TGfl ~ 10 — 30 . 



the stretching regime, which can last up to ten orders of magnitude in time. In this 
regime, the energy transfer into loops is negligible; this can be seen in figure 5.2(b), 
where initially less than one loop is produced per Hubble time. 

As the string velocity increases, energy losses to loop production become more 
and more important and eventually the network 'switches' rather quickly to the 
Kibble regime. Here, the fraction of the energy density in long strings that is 
transferred into loops per Hubble time is in fact larger than unity! Of course 
there is nothing unphysical about this, since the string network also gains energy 
during the time interval in question due to stretching — most of which is immediately 
converted into loops. Note that in this regime, although the fraction of transferred 
energy is constant, the number of loops produced per Hubble volume per Hubble 
time decreases: this is because the correlation length is 'catching up' with the 
horizon, so larger and larger loops (relative to the horizon ) are produced, given our 
assumption that during the friction-dominated epoch a ~ 1. 

In the opposite regime for high initial densities, the initial velocity is much 
larger, so loop formation is important right from the start and the Kibble regime 
begins immediately. Note that the first few orders of magnitude in time after the 
formation of an electroweak string network are the relevant period for baryogenesis 
mechanisms; we therefore believe that these results can shed some new light in this 
area. 

As we approach t eq , the network scaling switches again to a different regime, 
though now rather more slowly. This would be the matter-era analogous of the 
Kibble regime (L oc t 3//2 , v oc t 1 / 2 ), except that it is not particularly distinct given 
that the friction-dominated dynamics ends and the network evolves towards the 
final linear scaling regime. As we have already pointed out, the long-string network 
reaches the linear scaling regime about today. Nevertheless, it is still building 
up small-scale structure and, if our ansatz is valid, it will keep on doing so for 
another 20 orders of magnitude in time if gravitational radiation were the dominant 
decay mechanism (electromagnetic losses would intervene first for superconducting 
strings). 

It should be noted that although string dynamics is friction-dominated until 
t+, the same is not true for the corresponding term in the evolution equation for 
L (2.20) (this is of course due to the v 2 dependence of that term). In fact the 
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Figure 5.3: The approach to scaling for an axionic long-string network. Plots successively 
represent the exponents of the power- law dependence of L and V, V itself, and the ratio of the 
long string and background densities. The horizontal axis is labeled in orders of magnitude 
in time from the moment of string formation, and (c) and (d) are log plots; initial conditions 
are specified in table 5.2. 

expansion term in (2.20) is always dominant, although in the Kibble regime the 
friction and loop production terms are of the same order of magnitude. Also note 
that, as we said before, the loop formation term is proportionally more important 
in the Kibble regime than in the 'free' scaling regime — in fact, the (linear-log) plot 
of the ratio of the friction and expansion terms in (2.20) is indistinguishable (apart 
from the scale of the y axis) from figure 5.1(a)! 

Axionic and GUT strings reach the final scaling regime well before t eq , so we 
will separately consider the approach to scaling in each case and then the transition 
between the linear scaling regimes in the radiation and matter eras. 

Figures 5.3-4 depict the approach to scaling of an axionic string network. The 
discussion is analogous to the one for the early evolutionary stages of electroweak 
strings, except for the different orders of magnitude involved. The other difference 
is of course that the additional logarithmic dependencies affect all the transient 
scaling laws with the exception of L oc t 1 / 2 . Note that in all cases the effect of 
the logarithms is to increase the exponent of the power-law dependencies. The 
stretching regime can now last up to five orders of magnitude in time. 

C. Global and gauge GUT strings 

The approach to scaling for a global GUT string network is shown in figure 
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Figure 5.4: Loop production during the approach to scaling of an axionic string network. 
The horizontal axis is labeled in orders of magnitude in time from the moment of string 
formation, and plots represent the fraction of the energy in long strings converted into loops 
per Hubble time (a) and the number of loops created per Hubble volume per Hubble time 
(b). We have taken a ~ q (see (2.28)). 

5.5; as we mentioned previously, in this case loop reconnections onto the network 
are not important and hence the cases plotted are representative of all physically 
meaningful initial conditions. In this case the period of friction-dominated evolution 
is much shorter, and so the stretching and Kibble regimes, although still clearly 
distinguishable, are not very definite. Again, the effect of the logarithmic corrections 
for global strings is to increase the power-law dependencies. It takes a global GUT 
string network about six orders of magnitude in time to reach the scaling regime. 

Now, let us concentrate on the more interesting case of gauge GUT strings for 
which the approach to the linear regime is shown in figures 5.6-7. Here the epoch 
when friction dominates the dynamics lasts less than three orders of magnitude in 
time, and the linear regime is reached about six orders of magnitude after string 
formation. Hence the linear scaling regime is approached faster than previously 
estimated and, more importantly, strings become relativistic at a very early epoch. 
Consequently, energy losses to loops are relevant at all times (see figure 5.7). One 
should also notice the similarity between the plots of the exponent of the power-law 
dependence of the lengthscale L (labeled (a) in 5.6) and of the fraction of the energy 
density in the form of long strings converted into loops per Hubble time (labeled (a) 
in 5.7). This is evidence of the fact that loop formation is the crucial mechanism 
for cosmological string network evolution. It should be emphasized that with our 
ansatz for ct, ot sc ~ 10 ^ and T^G/j ~ 65 x 10 ^ it takes about 4 orders of magnitude 
in time for a to evolve from a = 1 to a = a sc , and during this time the physical 
size of the loops formed changes by less than a factor of 3. This is in agreement 
with both simulations [4] and previous analytical estimates [6]. 

Note that for the given parameters the energy density in string loops in the 
linear regime is in fact larger than that in the long strings. In this case, our analytical 
estimate was g es t ~ 2.1 in the radiation era, while numerically we find g nU m ~ 3.26 
(see figure 5.7(c)). Furthermore, while with our simplifying assumptions one finds 
that in the radiation-era scaling regime loops decay at a time tdecay ~ 5.18t / after 
formation (see figure 5.7(d)), the exact result is tdecay ~ 6.22t/. This is explained 
by the fact that in the analytic estimate we have assumed that loops always have 
V£ = 1/ sqrt2, while numerically we assume that their initial velocity is that of the 
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Figure 5.5: The approach to scaling of a global GUT string network. Plots successively 
represent the exponents of the power-law dependence of L (a) and V (b), V itself (c), and 
the ratio of the long string and background densities (d). The horizontal axis is labeled in 
orders of magnitude in time from the moment of string formation, and initial conditions are 
specified in table 5.2. 



long string network at the moment of formation. These discrepancies get smaller 
as we decrease the parameter a (hence the loop size at formation). 

It should also be pointed out that although less loops are produced in the 
friction-dominated regime, they live much longer — up to four orders of magnitude 
in time for loops produced just after the network forms. In fact, with our choice 
of parameters (note that a is close to the numerical upper limit) no loop decays 
during the friction-dominated epoch, so the energy density in loops relative to that 
in long strings grows to a value much larger that the final value in the linear regime. 
Hence if we had included the loops present at the moment of the network formation 
(which account for about 20% of the string length) we would find that loop density 
is always a non- negligible fraction of the long-string density. As one would expect, 
the loop population takes a longer time to reach the scaling regime compared to 
the long-string network. 

Finally, we discuss the transition between the radiation and matter era linear 
scaling regimes (see figures 5.8-10). Using our ansatzes for k and a, and assuming 
a smooth change of c (although this is strictly not necessary given the relatively 
weak dependence of the scaling properties on c) we can match all results the from 
numerical simulations [4] and our own analytical estimates (see section 4c) with 
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Figure 5.6: The approach to scaling of a gauge GUT string network. Plots successively 
represent the exponents of the power-law dependence of L (a) and V (b), V itself (c), and 
the ratio of the long string and background densities (d). The horizontal axis is labeled in 
orders of magnitude in time from the moment of string formation, and initial conditions are 
specified in table 5.2. 

an error of less than 10%. This is a significant achievement given the fact that in 
previous analytic models the discrepancies were typically of order 100%. 

Note that this is a much slower and smoother process than previously esti- 
mated, extending for about eight orders of magnitude in time. This is the reason 
why Bennett & Bouchet [4], having studied this transition with a series of numeri- 
cal simulations covering a range of five orders of magnitude in time, do not see the 
scaling parameters "going flat". In fact, one can say that the string network never 
leaves the scaling regime — the exponent of the power-law dependence of L is never 
more than ten percent away from the scaling value. It is also interesting to observe 
that the string velocity increases just after t eq , before it decreases to the matter-era 
scaling value (see fig. 5.8(d)); this is due to the comparatively rapid change in the 
expansion rate. 

Note the effect of the logarithmic corrections for global strings in the evolution 
of the ratio of the long string and background energy densities (see figure 5.9). It 
should also be said that due to its specific parameter dependencies, the number of 
loops produced per Hubble volume per Hubble time is the quantity for which our 
temporary matching problem is more serious. 
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Figure 5.7: Evolution of the loop population for a gauge GUT string network. Plots 
successively represent the fraction of the energy in long strings converted into loops per 
Hubble time (a), the number of loops created per Hubble volume per Hubble time (b), the 
ratio of loop and long string densities (c) and the time of loop decay (in orders of magnitude 
after formation) for loops produced at each time (d). Note that (a,b,d) were obtained from 
our analytical estimates (see section 4C) while (c) is the (exact) numerical solution of (4.30). 
The horizontal axis is labeled in orders of magnitude in time from the moment of string 
formation, and we have taken a ~ 10 -3 and TG[i ~ 65 X 10 -6 . 

To conclude, it should be emphasized that this quantitative picture of GUT 
string evolution can lead to important modifications in the structure formation 
scenarios involving cosmic strings. 

6 Conclusions 

In this paper we have presented a detailed account of a recently-developed gen- 
eralized 'one-scale' model of string network evolution [8,16] where a 'characteristic 
lengthscale (the 'correlation length' or average inter-string distance in the case of 
long strings, the physical loop length in the case of loops) and the average velocity 
are the dynamical variables. 

The immediate benefit of this generalization is that one is thus able to properly 
describe string motion in friction-dominated contexts. As a consequence, this simple 
model is the first complete and fully quantitative study of the evolution of a string 
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Figure 5.8: Evolution of GUT long-string networks (both gauge and global) in the tran- 
sition between the radiation and matter linear scaling regimes. Plots respectively represent 
the exponent of the power-law dependence of L (a) and V (b), the ratio L/t (c) and ve- 
locity (d). The horizontal axis is labeled in terms of the logarithm of the scale factor (with 
d{t eq ) = 1); it spans the period between 10~ 10 t eq and 10 10 t eq . 



network and the corresponding loop population in both condensed matter (see [16]) 
and cosmological contexts. 

Notably, we established the validity of Kibble's scaling law for intermediate 
and light energy scale cosmic strings; furthermore, if the initial string density is 
sufficiently low, there can also be an earlier period during which the strings are 
conformally stretched. However, despite the relative growth in the string energy 
Pool Pb in the damped epoch, we showed that strings can never dominate the energy 
density of the universe. The model also predicts that electroweak strings are only 
approaching the linear scaling regime today, whereas GUT strings reach this scaling 
regime faster than previously estimated. Finally, we have shown that the transition 
between the radiation and matter era linear scaling regimes is a very slow process. 

Application of the model to string loops has also led to the determination of the 
evolution of the density and other relevant properties of the loop population. We 
can use this model to determine the loop density and other relevant properties at all 
times. We have found that for typical values of the parameters characterizing loop 
production and decay there is more energy density in loops than in long strings. 

These results can significantly affect some cosmological string scenarios. The 
insight gained on the friction-dominated epoch of string evolution also calls for a 
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Figure 5.9: Evolution of the ratio of long string and background densities for global (a) 
and gauge (b) GUT string networks in the transition between the radiation and matter eras. 
The horizontal axis is labeled in terms of the logarithm of the scale factor (with Cl(t eq ) = 1); 
it spans the period between 10 — 10 t e „ and 10 10 t e „. 




Figure 5.10: Evolution of loop characteristics for GUT string networks in the transition 
between the radiation and matter eras. Plots represent the fraction of the energy density in 
long string converted into loops per Hubble time (a) and the number of loops produced per 
Hubble volume per Hubble time (b). The horizontal axis is labeled in terms of the logarithm 
of the scale factor (with a(t eq ) = 1); it spans the period between 10 — 10 t eq and 10 10 £ e q; 
we have taken a ~ 10~ 3 and TG[1 ~ 65 X 10~ 6 . 

re-examination of the vorton problem [25] and of string-induced baryogenesis [26]. 
It should also be possible to obtain more accurate estimates of the contribution of 
axions to dark matte, and, of course, our results on the radiation-matter transition 
for GUT strings can affect structure formation scenarios. We hope to tackle some 
of these problems in the very near future. 

We have also pointed out that the study of the evolution of cosmic string 
networks has been an effort involving both analytic and numerical work. Our results 
concerning the friction-dominated epoch of string evolution make it clear that a 
numerical study of this epoch is long overdue, not least because (as was pointed 
out by Vilenkin [18]) including the frictional force is simplicity itself: it amounts 
to a redefinition of the Hubble parameter, as can be seen from (2.9,2.10). There is 
also a need for a more careful study of the loop population. We have made several 
predictions regarding the scaling ratios of the loop and long string densities as a 
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function of parameters characterizing the loop production and decay; these should 
be testable by (slightly modified) existing numerical simulations. 

Apart from the effect of reconnections back onto the long string network, which 
is non-negligible in some cases, the outstanding remaining problem is that of the 
small-scale structure seen in the simulations [4]. We believe that the parameter k, 
defined in (2.40), will ultimately provide a phenomenological means of introducing 
small-scale structure effects in this model. Some steps in that direction have already 
been taken in this essay — in particular, a possible solution for the famous 'matching' 
problem of Kibble's one-scale model has already become apparent. In particular, it 
seems likely that we can relate k to the effective ('renormalized') string energy per 
unit length p,, which is numerically seen to provide a good description of small-scale 
structure. These important issues are presently being considered. 
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